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ABSTRACT 



As an alternative to the covariant Ostrogradski method, we show that higher- 
derivative relativistic Lagrangian field theories can be reduced to second differential- 
order by writing them directly as covariant two-derivative theories involving Lagrange 
multipliers and new fields. Despite the intrinsic non-covariance of the Dirac's proce- 
dure used to deal with the constraints, the explicit Lorentz invariance is recovered 
at the end. We develop this new setting on the grounds of a simple scalar model 
and then its applications to generalized electrodynamics and higher-derivative grav- 
ity are worked out. For a wide class of field theories this method is better suited than 
Ostrogradski's for a generalization to 2n-derivative theories. 
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1. Introduction 

Higher-derivative (HD) theories have an old tradition in physics. Relativistic 
field theories with higher order Lagrangians historically range from Higgs model reg- 
ularizations [1] to generalized electrodynamics [2] [3] and HD gravity [4], and arise 
as effective low energy theories of the string [5] or induced by quantum fields in a 
curved background [6] . A procedure was later devised to reduce them, by a Legendre 
transformation, to equivalent lower-derivative (LD) second-order theories [7] where a 
subsequent diagonalization explicitly displays the particle degrees of freedom [3] [8] [9] . 

The validity of the formal Lorentz covariant order-reducing method adopted there 
has been checked in an example of scalar HD theories by a rigorous study of the 
phase-space [9]. In this procedure, a generalization of the Ostrogradski formalism to 
continuous relativistic bosonic systems (2n-derivative because of Lorentz invariance 
in the most usual cases) is carried out. In it, some of the field derivatives and the 
generalized conjugate momenta become, after a suitable diagonalization, new field 
coordinates describing the degrees of freedom (DOF) which were already identified in 
the particle propagators arising in the algebraic decomposition of the HD propagator. 

By using Lagrange multipliers, an alternative to the Ostrogradski method for me- 
chanical discrete systems has been proposed which allows to show the quantum (Path 
Integral) equivalence between the modified action principle (first order Helmholtz La- 
grangian) and the starting HD theory [10]. For relativistic field systems, we show that 
a similarly inspired procedure can be followed in which the multipliers let to write 
the HD theory from the outset as a second order (constrained) covariant one which 
lends itself to a particle interpretation after diagonalization. 

In this paper we implement this new setting by means of the use of Lagrange 
multipliers in a Lorentz invariant formulation of relativistic scalar (and subsequent 
applications to vector and tensor) field theories. In the following we will generically 
refer to them as HD theories (HDT). The Dirac method [11] prescribes the identi- 
fication of the primary constraints arising in the definition of the momenta. These 
constraints are added to the starting Hamiltonian by means of new multipliers, and 
then they are required to be conserved by the time evolution driven by this enlarged 
"total Hamiltonian" through the Poisson Bracket (PB). This may give rise to sec- 
ondary constraints, the conservation of which can in turn generate more secondary 
constraints. The process stops when the equation obtained is not a constraint but 
an equation allowing the determination of a multiplier. We then use the remaining 
constraints to eliminate the remaining multipliers and the momenta, ending up with a 
two-derivative theory depending on just its true DOF. Since the latter appear mixed, 
a diagonalization works finally out the independently propagating DOF. 

As long as time evolution is analyzed, the true mechanical Hamiltonian (i.e. the 
energy) of the system must be used. Then one cannot benefit of the compactness of 
the Lorentz invariant procedures introduced in [9], so we are initially forced to deal 
with non covariant 3+1 objects and face the diagonalization of larger matrices. The 
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relativistic invariance of the system becomes explicit only at the end of the process. 
Alternatively this invariance may also be maintained at an earlier stage by merging 
the constraints into covariant equations: in fact it is possible to find all the previously 
known 2nd class (non-covariant) Dirac's constraints by working out the whole set of 
covariant field equations (for fields and multipliers). 

Prom the methodological point of view, the new treatment of HD boson theories 
that we present here provides a sharp departure from the more traditional Ostro- 
gradski approach. Moreover, it is implementable and may prove advantageous in 
generalized electrodynamics, HD Yang-Mills and linearized HD gravity as well. On 
the other hand it lends itself to generalization to arbitrary n better than the Ostro- 
gradski method does. 

Our work focuses on the propagating DOF of these theories. Therefore we mainly 
consider the free part of the corresponding (HD and their equivalent LD) Lagrangians. 
Self interactions (derivative or not) and interactions with other (external) fields are 
embodied in a source term, namely a source j linearly coupled to the fundamental 
field. The source always contains a coupling constant. Thus for our purposes one may 
adopt two attitudes: either j is assumed to contain only (spectator) external fields or 
the whole source term is treated as a perturbation (and ignored at the zeroth order in 
the coupling we are interested in) if self couplings are present. In any case we retain 
an external source as a guide to find the suitable diagonalizing matrix. 

HD (free) fermion theories will always be (2n + l)-derivative to maintain Lorentz 
invariance. The constraint method presented here (as well as Ostrogradski's) can 
be adapted to encompass these theories, but the task is beyond the purposes of this 
paper. 

Again resorting to the scalar example, in Section 2 we review the general situation 
of HD theories in which the "square masses" may be the complex and/or degenerate 
roots of a polynomial, then justifying the custom of focusing only on theories with real 
non-degenerate squared masses. In Section 3 we treat n = 3, the case n = 2 being 
too much trivial for our illustrative purposes. The results regarding the extension 
to arbitrary n are presented in Section 4. Through the Dirac's canonical method 
followed, they guarantee that all the constraints present in the theory have been 
properly taken into account. This validates the Lorentz covariant phase-space method 
described in Section 5 which we adopt in the following as the suitable shortcut for 
practical uses. Section 6 discusses, with the appropriate mathematical formalism, 
the application of our constraint method to HD (Podolsky) vector field theories and 
generalizations together the most interesting case of HD (linearized) gravity. We 
summarize our results in the Conclusions in Section 7. 

In Appendix A we show the pure algebraic character (i.e. absence of space 
derivatives) of the secondary constraints. Also for the general 2n-derivative case, 
in Appendices B and C we display the steps leading to the diagonalization of the 
DOF. In Appendix D, when squares of the multipliers occur, we show the equivalence 
between this seemingly constraint method and Ostrogradki's. 

Throughout the paper the Minkowski metric (H ) is adopted. 
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2. HD Scalar Theories 

The theories of arbitrary (finite) differential order have a distinctive spectrum 
with respect to the usual second order theories describing physical systems. It gener- 
ally includes peculiar ghostly states and non-particle solutions. Although this feature 
is not exclusive of HDT's (one may devise second order theories with this sort of com- 
plications), it is almost inherent to them. The exception corresponds to very special 
cases like HD gravity made up only with powers of the scalar curvature. There is no 
escape however for simpler fields. 

The simplest HDT's are linear theories of scalar fields 

N c 

C =-^q N { U )^-j^ , (2.1) 

where c N is a constant with suitable dimensions and q N (□) is a monic real polynomial 
in □. A detailed study of the theories (2.1) in which q N has only real non-degenerate 
roots can be found in [9]. In these cases, the theory (2.1) is equivalent to a LD one, 
made up of the alternate sum of N Klein-Gordon (KG) free Lagrangians where the 
sum of the respective real fields couples to the source j. This states that all the DOF 
are particle-like and couple to the same source. 

However this is a privileged case of (2.1) as long as q N may generally have de- 
generate and/or complex roots 

qJP) =(□ + mi) Rl •••(□ + m 2 r ) Rr (n - M 2 ) Tl •••(□- M 2 ) Tt 

[(□ - fii)(n - cii)} Cl •••[(□- n c )(n - o c )] Cc 

where the correspond to physical masses, Mj are tachionic and fij are complex, 
Ri, Ti and Cj being their respective degeneracies. 

Non degenerate complex roots can be handled formally with the methods devel- 
oped in [9]. Following them, it can be seen that 

c 4 = -j<P(n-mn-n)<p- 3 <p 

is equivalent to 

c 2 = ^(n-n)[ v ( n -n) v -p( n -ti)<p]-j( v + <p) , (2.2) 

where ip = ipi +if2 is a complex scalar field. However this is not a real diagonalization 
because (2.2) cannot be diagonalised in terms of real (or even complex) scalar fields 
built linearly out of <pi and </? 2 and with real square masses. 

Also the degeneracy makes it impossible to interpret the roots as the square 
masses of free particle-like (even tachionic) states. In fact, already the simple case 
</>(□+ m 2 ) 2 <p cannot be reduced to a LD theory of free real (or complex) fields because 
the HD propagator l/(d + m 2 ) 2 is not a linear combination of LD free propagators 
1/(D± m 2 ) . Otherwise stated, there are solutions to (□ + m 2 ) 2 (j) = that are not 
solutions to (n + m 2 )(p = 0, and hence cannot be expressed as superpositions of plane 
waves. 

Summarizing, only in the case of real non-degenerate square masses can the 
HD propagator be algebraically expanded as a sum of particle-like propagators, i.e. 
in terms of propagators displaying physical or tachionic masses (sign of the square 
mass), or even a non-physical norm (sign of the propagator). As customary in the 
literature, in the following we will limit ourselves to the study of these cases. 
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3. n = 3 theory 

We consider the six-derivative Lagrangian, 

£ 6 = -^0[lJ[2j[3j0-.7 , (3.1) 

where fx is an arbitrary mass parameter, [i] = (□ + mf) are KG operators, M = 
(12) (13) (23) , (zj) = m 2 -m 2 > for i < j , and mass dimensions [//] = 1, [M] = 6, 
[(f)] = 1, [j] = 3. We have overlooked the case n = 2 because it is too trivial for 
introductory purposes. 

As discussed in [9], (3.1) displays the general form of the free part of a higher- 
derivative scalar theory with non-degenerate masses mi , m 2 , ffi3 , the source term 
embodying the remaining self-interactions and the couplings to other fields. There we 
showed by a covariant Legendre order-reduction procedure that (3.1) is equivalent to 
the second-order Lagrangian 



^^-l^M^i + l^Un^-l^U^-J^i+^ + h) , (3-2) 

in agreement with what is expected from the algebraic expansion of the HD propagator 
in particle poles, namely 

l^~ 2 M ^~ 2 (23) | ^ 2 (13) ^ 2 (12) _ 



[1][2][3] [1] 

The physical meaning is that whenever a source j emits a field driven by a HD 
dynamics, it actually emits a linear superposition = 0i + 02 + 03 of LD (particle) 
fields. 

However £ 6 can be recast directly as a second-order theory with constraints, 
namely 

£6 = \m Hal* + - » 2 ) + A 2 (^ 2 - [3]^ 3 )H^3 , (3.4) 

where V> 3 = and Ai , A 2 are Lagrange multipliers, so that £ 6 depends on five fields. 
It is immediate to check that the equations of motion for Ai and A 2 yield ipi and ^2 
as (derivative) functions of $3. By substituting them in (3.4) and identifying -0 3 and 
one recovers the HD equation (3.1). Dropping a total time-derivative, in compact 
matrix notation, (3.4) reads 

C 6 = lq T JCq + ^ T M^- J T ^> , (3.5) 
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where the vectors \1/ and J, with components ipi , Ji , and the matrices /C and .M are 



1 



so that [V'i] = 1 i = 1, • • • , 5 ; Ji = jS; 
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(3.6) 
\ 
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/ 



Ai is an operator with space derivatives present in Mf 
The canonical conjugate momenta are defined as 



mi 



A. 



TV,. 



PC 6 

dipi 



(3.7) 



They are the components of a 5-vector IT for which one has 

n = /C^ . 



(3.8) 



Since K, is not invertible, not all the velocities ipi can be expressed in terms of 
the momenta and a primary constraint occurs, namely 



Ol = 7T5 — 7I~l = , 



(3.9) 



as consequence of Tr§ = 2M^3 = tti • There is only one such a constraint since the 
submatrix fC a b = 2M^ a b (a, b = 1, ... ,4) is regular. In the following, indices a, b, ... 

go from 1 to 4, while i, j, ... go from 1 to 5. The velocity ip$ is not worked out, and 
from (3.8) we have 



2M 



2M 



(3.10) 



and therefore 



2M / 

(I 



(3.11) 
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The Hamiltonian is 



1 • • u 6 • • 1 

U = n a ip a + 7T 5 ^5 - -ll)a^ab^b ~ ^^3^5 ~ -^iMijtfjj + jlf) 3 , (3.12) 

where ip a must be substituted according to (3.11). Then the dependence on ip 5 cancels 
out and we have 

1 2M 1 

H = -^—^-KaKabKb ~ ^i-Mijlfjj + Jilpi . (3.13) 

Since not all of the five momenta 7Tj are independent because of the primary 
constraint (3.9), only four momenta appear in (3.13) together with the five fields ipi . 
The "total Hamiltonian" , with five independent momenta, accounting for this is 

H T = H + (n 1 , (3.14) 

where £ is a Lagrange multiplier. 
The stability of Q\ requires 



tt 1 = {n 1 ,n T } PB = n 2 = ^((13)iJs-^4 + ^2)=0 . (3.15) 
This secondary constraint yields 

^4=^-^3+1^2 ■ (3.16) 

Further secondary constraints stem from the ensuing stability conditions 

n 2 = {n 2 ,n T } PB =n 3 = (13)^-^2 + ^4 = , (3.17) 

so that 

Tr^^-Myn , (3.18) 

and again 

n 3 = {n 3 ,n T } PB = n 4 = 

h ■(-(13)(23)^-^ 2 (13)^ 2 -^V 1 + ^V 5 ) =0 , 
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(3.19) 



2M 



(once (3.16) has been used), from which one gets 



+ . (3.20) 



The next constraint, after using (3.18), gives 



O 4 = {O4,HT} PB =O5 = (13)(12)7r 1 -^ 2 (13)7r 2 -/7r 3 + 2^:^ 4 C = , (3.21) 

and Q can be obtained as a function of 7i"i , i\i , and , thus bringing the generation 
of secondary constraints to and end. 

1~Lt being quadratic in t/j's and 7r's, guarantees an alternance of linear constraints 
involving the fields and the momenta. In spite of the occurrence of space derivatives 
in M., they cancel out and the constraints are algebraic. From this set of constraints, 
the multipliers 1(^4 and ^5 , together with their conjugate momenta 114 and tt^ , can be 
worked out. 

The Lagrangian (3.5) can be expressed in terms of the independent variables 
ip a (a = 1, 2, 3) . Notice that implementing these algebraic constraints in C 6 does 
not modify the second differential order already attained in (3.5). One obtains 

£ 6 = ip a fc a pij>0 + il> a M a pil>p - j^z , (3.22) 

where 



= ^(X*/3 + K-olbNbP + NaAK-Ap) = ^ ( 1 ,.V 2 '_, I ■ (3.23) 



( f -Ml \ 






1 


1 


(13) 




M 2 


(13) 


(13) (23) 


M 2 


M 4 



2M 



I? (13) -Ml -^M? 



V-M| -<i|>Mf -^P-MIJ 



(3.24) 



with a, /?, . . . = 1, 2, 3 ; A, B, . . . = 4, 5 ; and 



1 < 13) 



Na/3 — I ^ (13) (13^(23) I (3.25) 
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that allows to embody (3.16) and (3.20) in the closed form 



i>A=M A ^p ■ (3.26) 

The symmetric matrices IC and M. can be simultaneously diagonalized by the 
regular transformation 



= , (3.27) 



where 




K aP = _ Uf> o • (3-28) 



The reason for this property will be made clear later on. 

The 3rd. row of the non-orthogonal matrix 1Z in (3.28), has been chosen so as 
to yield the source term in (3.2), in which the source couples to the sum of the LD 
effective fields. The remaining six elements are uniquely determined by requiring 1Z 
to diagonalize K, and M.. 

The diagonalized matrices are 



^Mn=^ diag (-M ? ^,Ml™-Ml^p) . (3.30) 

so that (3.22) finally writes as C 2 in (3.2). 

This shows again the Ostrogradski-based result [9] of the equivalence between the 
six-derivative theory (3.1) and the LD version (3.2) that reproduces the propagator 
structure (3.3). 



9 



4. Theories with arbitrary n 

The general Lagrangian 

£ 2n = ~^[l][2]...[n]0-j* , (4.1) 
where M = Y[ (u) ? an d <i = n(n — 3) + 2 for dimensional convenience, can be dealt 

i<j 

with along similar lines. The 2-derivative constrained recasting of (4.1) is 

C 2n = ^[-^»[l]^i + Ai(^i-[2]^) + ... + A n _i(^ B _i-[n]^ n )]-^„ , (4.2) 

with ip n = (f) , and Ai, . . . , A n _i being Lagrange's multipliers. In order to have a more 
compact notation we define 



ipa = V 2(n a) ipa ; a = l,...,n 
ipA = A t_2a Aa ; A = n + a;a = l,...,n — 1 
so that [V>i] = 1 (i = l,...,2n-l). Then 



(4.3) 



— ■ — ■ 

with Ji = j5i n , and the (2n — 1) x (2n — 1) matrices /C and M. are given by 



(4.5) 



l . n ( n — i) 

with a = - 2M — . This "mass" matrix contains again space derivatives. Here and in 
the following an underlined index means that Einstein summation convention does 
not apply. The canonical conjugate momenta are now 

dC 2n 



i.e., in closed notation, 



Defining the matrix K, 



II = /C* . (4.7) 



1 



IC ab = -IC ab (a,6=l,...,2n-2) , (4.8) 



one sees that det /C 7^ 0, while det /C = 0. This means that we only have one primary 
constraint, namely 

fix EE 7T 2n _i - 7Ti = . (4.9) 
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Then ip2n-i is not worked out, while ip a (a = 1, . . . , 2n — 2) can be expressed in 
terms of n a and ip2n-i- The first 2n — 2 components of eq.(4.7), namely 



n a = aK, ab ip b + aSan4>2n-l (4.10) 

give 

1p a = -K-ab^b ~ S a i1p2n-1 ■ (4.11) 
(7 

After checking that the terms in ip2n-i cancel out, the Hamiltonian has the simple 
expression 

1 ' 1 

ft = Tf^aKabKb - ^tpiMijipj + jip n . (4.12) 

In Ti only 2n — 2 momenta n a occur against 2n — 1 fields ipi , because of the 
primary constraint (4.9). One may restore the dependence on 2n — 1 momenta by 
introducing the "total Hamiltonian" 

H T = H + (ni , (4.13) 

where ( is a Lagrange multiplier. 

From the stability condition on Hi , a cascade of secondary constraints follows, 
eventually ending with an equation that determines the value of (. We outline here 
the steps closely following the lines of section 3. 



(In) 

n 1 = {n 1 ,n T } PB = n2 = o ^ n+1 = Vn-i + ^rVv* • (4.14) 

Then 



2 = {0 2 ,^ T } pB = 3 = /U 2 7r 2n _ 2 + (ln)7r 1 -^ 2 7r 2 (l-5 n2 )-2(7C5 n2 = . (4.15) 

If n = 2, eq.(4.15) gives £ in terms of t\\ and 7r 2 , and the cascade stops here, but if 
n > 2 it yields 

(In) 

7T2n-2 = — W-TTl + . (4.16) 

The next step is ^3 = {O3, 7~Lt} PB = ^4 = 0, which together with (4.14) gives 



V>n+2 =Vn-2 + ^ 2 -((l^) + (2,n-l))^ n _ 1 + ^(ln)(2n)V'n , (4.17) 
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and, proceeding further, we obtain for the momenta 



4 = {^4,'Ht} pb =^5 = /U 4 7T2n-3 ~ (in) (1,71 - 1)tTi + 

+ fi 2 ((ln} + (2,n-l})n 2 -fi 4 n 3 (l-5 n3 )-2a^C5 n3 = , 

where (4.16) has been taken into account. Again, if n = 3, the process stops here and 
we have reproduced the results of section 3. If n > 3, eq.(4.18) yields 

7T 2n _ 3 = ^(ln)(l,n-l)7r 1 -^((ln)(2,n-l))7r 2 + 7r 3 , (4.19) 

and the process goes on. 

For illustrative purposes, we complete here the steps that cover the case n = 4. 
fi 5 = {Q 57 n T } PB = ^6 = 0, yields 



tfj n+3 = ^ n _ 3 + -J «3, n - 2) + (2,n - 1) + (ln»^ n _ 2 + 

+ -1 ((2, n - 1)(3, n - 1) + (ln)((2, n - 1) + (3n)))^„-i+ (4.20) 
+ 4(l™>(2n)(3n)Vn , 

and 06 = {^6, T~Ct} pb = ^7 = gives £ in terms of 7Ti , tt 2 , 7r3 , and 114 . 

In general, for a fixed n, the quadratic dependence of H on 7Tj and ^ , together 
with the primary constraint fii , leads to a set of secondary constraints fifc that splits 
in two classes according to k being even or odd. A constraint fi 2 j (j = 1, . . . , n — 1) 
is a linear combination of tpi and gives ift n +j i n terms of ift n , . . . ,ij) n -j. A constraint 
2j _i (j = 2, . . . , n — 1) is linear in m and gives 7r 2n _j in terms of 7Ti, . . . , 7Tj . Finally, 
2n _i fixes the value of ( and stops the process. 

One can prove that the constraints on the momenta 2n _i do not contain space 
derivatives, even though the elements of M. involved in their calculation contain the 
Laplacian operator. This will be shown in the Appendix A. 

Like in (3.26), we take 

V>a = M A ^a (4.21) 
with indices a, (3, . . . = 1, . . . , n and A, B, . . . = n + 1, . . . , 2n — 1. 
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Now Af is a (n — 1) x n numerical matrix whose three first rows can be read 
from (4.14), (4.17) and (4.20). In order to write the elements of A/", it is useful to 
introduce the indices s = A — n, (s = 1, ...,n — 1) , that labels the rows of A/", and 
p = A + (3 — 2n = (3 + s — n, (p = 1 + s — n, s), that indicates the odd-diagonals of 
j\f. Then, with P^+j 1 = _/V n+s , n _ s+p , one has: 



Pg} = for p<0 ; P? +1 = l ; P| +1 = -I £ (ji, n + 1 - h) 

^ s-1 s-j'i 



^ Jl=l J2 = l 



^ s-p+1 s-p+2-ji 

P pTi=^ E^' n+1 -^ (ji+h,n + 2-j 1 -j 2 ) 

^ 31=1 J 2 = l 

s-p+l-j 1 -...-ji- 1 

Z O'l + -■- + Jz, n + Z — ji -...-ji)... 

31 = 1 

S + l-j'l-...-j'p-2 

Z 0'i + - + i - ji - ••• - Jp-i) 

3p-i=1 
s-ji----j p -i 

Z 0'i + - + ip» 7i + p-ii---jp) 

jp=l 



1 2 3-ji 

= "2(^1) E^' n+1 - ^> Z O'l + J2, n + 2 - J! - j 2 ) . . 

^ 31=1 J2 = l 

s-l-ji- — -3s-3 

Z O'l + ••• + 3s-2,n + s - 2 - ji - ... - j s _ 2 ) 

Js-2 = 1 
S-jl-...-j B -2 

Z O'l + ••• + j s -i,n + s — 1 — Ji — ••• -ja-l) 

js-l=l 

P s s + 1 = -^(ln)(2n)...( S n) 



and (4.21) can be rewritten as 



(4.22) 



i; n+s = J2Ppt^n-s +P (4.23) 

p=0 
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The proof of (4.22) and (4.23) is given in Appendix B, were the relation 



P 8 p tl = P pT S for p>0 (4.24) 

is also proven. Taking P a p = Pp+i , with a = s + 1 , (3 = (p+1) — (s + 1) +n = p — s+n, 
and Pi n = 1 , P\p = for (3 < n, we can consider P a p as a n x n matrix, that is 
nothing more than the matrix M enlarged with a first row, which is symmetric as a 
consequence of (4.24). 

After using the constraints (4.23) to keep only the independent variables, the 
Lagrangian again is 



-2n 



Ipafcaplpp + 1p a M a plp{3 ~ jlpn ■ (4.25) 



The n x n matrices /C and M. have the same structure in terms of K, , M. and M 
given in (3.23) and (3.24). Taking into account (4.22), one checks that 

^n(n-l) ^n(n-l) 
K/CtP = 2M ^-n+a = 2M P a/3 = ^ Q (4.26) 

is the generalization of (3.23). For M. one has two contributions, with and without 
the operators Mf . The first one can be written, for a < f3 as 

n(n-l) 

M aP = - !t——MlP a p , (4.27) 

which displays for f3 < a. 

The contribution without is given by 

M a f3 = if a + [3 < n 

„n(n — 1) 

MaP = 2M~ ^ P ^ +1 if a + /3>n; a,[3^n (4-28) 
■Map = if a or (3 = n . 

The formulas (4.27) and (4.28), are the generalization of (3.24) for arbitrary n . 

The diagonalization of (4.25) will be accomplished, as in (3.27), by a n x n real 
matrix 71. We again impose 1Z n [3 = 1 , {(3 = l,...,n) to ensure that the current 
couples to each one of the degrees of freedom. The requirement of simultaneously 
diagonalizing K, and A4, yields n(n— 1) quadratic equations that determine the n(n— 1) 
remaining elements of TZ. The existence of such a regular 1Z with real elements is 
guaranteed by the underlying Lorentz covariance. In fact, the constraints used to get 
from (4.4) to (4.25) are the (covariant) ones on the fields. In the next section we 
will see that they can be obtained also from the (covariant) field equations of (4.2). 
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They can be implemented on this Lagrangian which can be then diagonalized directly 
avoiding the decomposition 3+1. The diagonalizing matrix turns out to be 



K a( 3 = 1 ; (a = n) 

n aP = (-iy-«Li- 2{n - a) (P,a + l}(P,a + 2}...(P,n} ; (J3 < a < n) (4.29) 

n a(3 = o ■ («</?), 

as is shown in the Appendix C. Of course, for n = 3, this 1Z is just (3.28). 

The diagonalized matrices are now the generalization of (3.29) and (3.30) 



2M V u2(n-l) 'V I u2{n-l) 

(4.30) 

_, (12)(13)...(l,z)(z,z + l)-(i,n) (ln)(2n)...(n-l,n ) x 

' 1 J ^(n-l) '-' fj2(n-l) 



^2(n-l) 



, _, +lM2 (12)...(lz)(z,z + l)...(zn) M2 (ln)...(n-l,n) 
'"'^ 2 * / u 2 ( n ~ 1 ) n / u 2 ( n ~ 1 ) 

With = TZ a/ 3(f)f3 , and integrating by parts, the final Lagrangian reads 



(4.31) 



£ 2n = (-l) n ^i Wi + (-l) n -^4r02[2]02 + • • • 

W W (4.32) 

• • • + (-l) n_i+1 ^^W^ + • • • " l^nlnUn - j{<h + • • • + 0n) , 

where (i) = —^-rjjM Y[ r^jyj- This is the result expected from the covariant Ostro- 
gradski method shown in [9]. 
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5. The covariant phase space 



The constraints on the fields and momenta obtained by the Dirac's method de- 
fine a subset of the phase space in which the solutions to the equations of motion are 
contained. In non pathological situations, there is a one to one correspondence be- 
tween points of this subspace and these solutions [12]. In the system under study we 
have found algebraic constraints on the fields that can be derived also directly from 
the equations of motion. They define subsets of the configuration space in which the 
trajectories of the system lie. Now we easily see how these constraints are derived 
(we put j = for simplicity). 

From (4.4) the following equations of motion are obtained respectively for the 
fields ip n +i 5 ipi (1 < i < n) and ip n . 

fj, 2 ipi = {i+ ljipi+i 

/J 2 1p n+i = [z]V>n+i-l (5.1) 
= [l]Vl + [*2n-l 

From them it is possible to work out the n — 1 fields ipi (n < i < 2n) in terms 
of the n independent fields ipi (1 < i < n) . This is done by using the first of the 
equations (5.1) to obtain 

Dipi+i = H 2 ipi + m 2 +1 ip i+1 (1 < i < n) (5.2) 

Then iterative use can be made of the second of the equations (5.1), taking (5.2) 
into account in each step, to obtain the algebraic constraints on the fields (4.14), 
(4.17) and their general form (4.23). The action (4.25), already depending only on 
the independent fields, explicitly reads: 

1 n—l s 

£ 2n = g hmi]v>i + E E p;i^n-s +P (^s -[s+ i])v s+ i] (5.3) 

s = l p=0 

where the non-diagonal structure is apparent. 

The diagonalization matrix (4.29) is suggested by the equations of motion stem- 
ming from (5.3) and directly derivable also from (5.1): 



^""^i = P] • • • N 

^ 2(n - 2) V>2 = [3]---HW 



(5.4) 



\x 2 ip n -\ = {nj ip. 



n 
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Defining n new fields a (1 < a < n) , with ip n = 4> a and obeying [a]0 a = , the 

a=l 

equations of motion (5.4) become 

^ 2(n - 1} Vi = (2, 1) • • • (n, 1) 0! 

= (3, i) . . . ( n? i) ^ + (3, 2) • • • (n, 2) 02 
(5-5) 



H 2 ipn-i = (n, 1) 0i H h (n, n - 1) 

= 01 H h 0n 



From these equations it is immediate to read out the elements of the diagonalizing 
matrix (4.29). 

The procedure implemented in this section provides a useful shortcut to get the 
diagonalized theory. However it is worth noticing that one can not avoid the arduous 
path followed in section 4 because it guarantees the consistency and the stability of the 
constraints. Moreover, the Dirac's procedure (when it exists) yields the structure of 
the reduced phase space (null set of all the constraints, modulo gauge, equipped with 
the induced non-degenerate symplectic form). The knowledge of the reduced phase 
space is necessary, for example, when one is interested in using path integral methods 
because it is on this space where the canonical quantization takes place. Then the 
covariant method presented in this section complements the results obtained in section 
4 besides being specially suitable for the DOF diagonalization. 
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6. Applications to other theories 

The constraint method we have developed for scalar theories can be implemented 
for HD vector and tensor theories as well. However the canonical 3+1 procedure 
described in Section 4 is not convenient for practical applications, mainly when more 
complex theories as these are considered. The simplified covariant version introduced 
in Section 5 gets rid of the intermediate constraints 2 j+i on the (non Lorentz- 
covariant) momenta and obtains the useful ones Q,2j on the (covariant) fields by 
deriving them directly from the (always covariant) equations of motion for all the fields 
(original HD, auxiliary and multipliers). Moreover when gauge symmetries occur there 
are further constraints, namely the first class ones, that make more cumbersome the 
process of working out the Lagrange multipliers. In any case, the covariant phase- 
space is the method of choice. 

HD vectors and differential-forms 

In the case of HD vector theories one can incorporate gauge symmetries to the 
discussion. A typical example is the generalized QED proposed by F.Bopp and 
B.Podolsky [13] 

£ = ~\f^ - A- 2 F, vU F^-j^ . (6.1) 

The structure of its constraints, in the gauge-fixed case, has been studied in [3] by 
a canonical forcefully non-covariant analysis carried out on the Ostrogradski-based 
order-reduction procedure. 

A recasting of a higher-derivative gauge-invariant Yang-Mills theory as a two- 
derivative one by means of constraints has been done in a non-covariant 3+1 way 
[14]. Of course, in the non gauge-fixed theory, second class constraints arise which 
coexist harmlessly with the first class relevant ones. However we are interested in 
keeping the explicit Lorentz covariance and, to this end, we make use of the same 
algebraic manipulations introduced in section 5. 

The theory described by (6.1) can be considered a member of a wider family of 
theories which share the same general structure of the scalar models (4.1) 

£2?1 = ^M^ A * ,11|2, - W ^ JA4 ' (6 - 2) 

where (p is a s-form (1-form in the generalized electromagnetism) , A is the exterior 
product, * the Hodge dual (associated to the Minkowski metric), and now the [z] 
operators must be understood as combinations of the exterior derivative d and its 
dual S such as 

[i]=W + m?. 
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Notice that if m n = the Lagrangian (6.2) describes a gauge theory with gauge 
transformations i— > + dA , where A is an arbitrary (s — l)-form. 

With all these considerations, the field equations derived from (6.2) have the 
same look of equations (5.1) for the scalar model. Moreover, all the arguments that 
follow this equations in section 5 apply in the same way to the (6.2) HD theory that 
depends now on differential-form fields. Notice that, despite the gauge invariance of 
the theory, the order reduction method and the linear redefinition that diagonalize 
the DOF in the Lagrangian work exactly as in the scalar case. 

HD Gravity 

The covariant Ostrogradski order-reduction of the four-derivative gravity leads 
to a two-derivative equivalent in which the particle DOF can be fully diagonalized in 
both the Diff-invariant case [8] and with a gauge fixing [15]. 

The constraint technique for the order-reduction of a pure four-derivative confor- 
mally invariant gravitational Lagrangian has been already used in a 3+1 non-covariant 
form [16], where further first class constraints from Diff-invariance occur. In a covari- 
ant treatment and for the general case including also two-derivative terms [17], a 
seemingly similar method is adopted where in place of the Lagrange multiplier a less 
trivial auxiliary field featuring a squared (mass) term is used. A little work shows how- 
ever that this method is identical to the covariant Ostrogradski' s [9]. We illustrate 
this on the grounds of the scalar model in the Appendix D. 

The covariant constraint method introduced in this paper provides a new ap- 
proach. The most immediate application in higher-derivative gravity regards the 
linearized theory, usually considered when analyzing the DOF. Take for example the 
four-derivative Lagrangian 

C = y^g[aR + bR 2 + cR liV R t " / ] . (6.3) 

The linearization around the flat Minkowski metric, namely = rj^ + , 

simplifies it to 

£ = \K^ v ^r^ + r^Q^ a ?(b, c)r aP , (6.4) 
where r^ v comes from the linearization of the Ricci tensor 

ffxu = T^-^lu aph ^ = — [fl h^vol ~l~ h u,^,a Oh^^ tl a,/j,v] i 

and G^ v a/3 and Q^ v a/3 are numerical matrices: 

qhvoc/3 — \ ^q^ u fj a ^ — rf a rf^ — rj^^rj^ 01 ^ 
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The field equation for h^ u can be written (omitting indices) in terms of the objects 
above in the form 



GTZ + -GQ^G] 



GKh = , 



(6.5) 



Its straightforward to show that the general solution to (6.5) is 



h = hi + h 2 , 



(6.6) 



where hi and h 2 satisfy, respectively 




(6.7) 



Omitting indices, the order-reduction of the theory by means of a Lagrange 
multiplier yields the two-derivative local Lagrangian 



where hi^ u is a new field and h^v is the multiplier. Of course, because of the 
Diff-invariance, first class constraints will remain when the Dirac procedure is carried 
out. However, it is possible to avoid the use of Dirac method if one looks directly at 
the space of solutions of the field equations. The field equations for the Lagrangian 
(6.8) are 



The last equation leads to G1Zh 2 = [i 2 hi. Then, the second one expresses algebraically 
hs in terms of hi and hi'- 




(6.8) 




(6.9) 



h 3 



a 



^h 2 - GQGhi . 



2/j : 



Using this constraint in the Lagrangian (6.8) we obtain 



£ = 



a 



^h\Gnh 2 + 2h\GQUh 2 - fi 2 h\GQGhi . 



(6.10) 



2/j' 
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The last step is again a diagonalization. The hint is given by (6.6)-(6.7), identifying 
h = h 2 . This leads us to introduce new fields hi, h 2 




GQ~ 1 Gh 2 



in terms of which the Lagrangian is 



2^ 2 



a 



hi Gfthi - 



2^ 2 



a 



h 2 GUh 2 - 




h 2 GQ- L Gh 2 , 



where now is clear that hi propagates the massless graviton DOF whereas h 2 carries 
a spin 2 with square mass —a/c and a scalar with square mass a/2(c + 36). 

7. Conclusions 

We have shown how to deal with 2n-derivative relativistic scalar theories by 
writing them directly as second-order constrained Lagrangians with more fields and 
suitable Lagrange multipliers. The corresponding canonical conjugate momenta are 
subject to primary constraints, whose conservation in time gives rise to a finite chain 
of secondary constraints according to the Dirac's procedure. Though expected, a non 
trivial result is that these constraints, later used to extract the final DOF, are purely 
algebraic relations that do not involve the space derivatives. 

Once the constraints have been implemented, we are left with a second-order 
Lagrangian for the DOF of the system. We have performed explicitly the diagonal- 
ization for n = 3, reproducing the result obtained in [9]. Then the procedure has 
been generalized to arbitrary n, namely (4.29). This explicit possibility constitutes a 
definite advantage over Ostrogradski's method. 

The applications to more interesting theories like HD generalized electrodynamics 
and HD Diff-invariant gravity illustrate also the fact that the order-reducing methods 
used in the literature fall in two categories: the one based in the covariant Ostrogradski 
and the one based in the constraints by Lagrange multipliers. The methods based on 
auxiliary fields with a quadratic term, which may look like a variant of the multipliers, 
actually belong to the first category and have no obvious manageable extension beyond 
the four-derivative order. 

In vector and tensor field theories where gauge symmetries occur, the correspond- 
ing first class constraints live together with the second class ones worked out in this 
paper and survive the order-reducing procedure as long as gauge fixings are not con- 
sidered. The method may then prove useful for a detailed analysis of the constraints 
from gauge (or Diff-)invariance in these HD theories, chiefly of the fate of the scalar 
and vector constraints of Hamiltonian gravity. 
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Appendix A 

We prove, by induction, that the constraints fi 2 j (j = 1, . . . , n — 1) involving 
the fields, do not contain space derivatives because the Laplacian operators cancel 
out. 

One first sees, by inspection, that this statement is true for fi 2 : the calculation 
leading to (4.14) is 

2 = a [/iV„-i + {Ml - Ml)ij) n - //Vn+i] = , 

(A.l) 

where the cancellation of the Laplacian operator is apparent, i.e. 

Ml-Ml = m\-ml = {\n) , 

(A.2) 

and obviously no summation is understood for repeated indices. Then, let us suppose 
that, after taking into account the preceding constraints, one has that in the constraint 

^2a = O [^Ipn-a +aiV>n-a+l +a 2 1p n -a+2 + - ■ . + a Q _i^ n _i +a a l(j n - fJ, 2a ip n+a ] = , 

(A.3) 

for a = 1, . . . , n — 2, the coefficients a\, . . . ,a a are real numbers, as are those found 
in (4.14), (4.17) and (4.20). We now prove that this is also true for 2a+1 . In fact 

^2a+l = /U 2a 7T2n-a-l + O-Y^ln-a, + Cl2^2n-a+l + ■ ■ ■ + aa-l7I~2n-2 + a a 7Ti — 

- ^ 2 7T a+ i(l - 5 n , a + l) - 2(T^ 2 C<5n,a+l , 

(A.4) 

from which 

^2(a+l) = ^2a+l = { ^2a+l , } pB = 

= a [-^ 2 °M 2 _ a V'n-a - fllM n 2 _ a+1 ^_ a+ l - ... - tta-iM^^n-l- 
- a a M 2 V>n + ^ 2Q M 2 +Q ^ n+Q + ^ 2 (° + 1 Vn-a-l + A*VV>n-a + ...+ 

+ n 2 a a -\^ n -2 + /j 2 a ol tjj n+1 - /i 2(Q+1 V n+a +i] = . 

(A.5) 
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The crucial point now is that, when working tp n + a out of (A. 3) and substituting 
it in (A. 5), only differences of squared masses Mf occur as in (A. 2), thus cancelling 
out the operators A. Then, by substituting also ip n +i from (A.l), one gets ip n+a+ i as 
a sum of linear terms in ip n7 ipn—i, . . . , ip n — a —i 1 the coefficient for the last one being 
the unity. This ends the inductive proof. 

Appendix B 

We proof (4.23) by induction. First note that for s = 1 it is nothing more than 
(4.14). Assume then that it holds for i/j n + s with s < n — 2. From the Lagrangian 
(4.4), one has that the equation of motion for ip s +i , s < n — 2, is 

H 2 ip n+S+1 = {s + l}ip n+s = (□ + m 2 s+1 )ip n+s 

(B.l) 

that, with (4.23), gives 

s s 

^tpn+s+l = Pp+lOlpn-s+p + m 2 s + i ^n-s+p , 
p=0 p=0 

(B.2) 

but the equation of motion for ^n-s+p-i is 



D^n-s+p — A* VVl— s+p— 1 ITln — s+plpn — s+p j 

(B.3) 

so 



^ +s+1 = ± {8+1 ' n 2 s+p) p;t^n- s+ p + ± pgfr 



n— s+p— 1 



p=0 

_ (s + l,n) ps+1 

,2 



s-1 



ptiHn + E 

p=0 



(S + 1,71 -S+p) pa+ i ps+ l 



V'n — s+p + 



P =l 



{s + l,n-s + p-l) +1 +1 ~ 



/7/ 



P+l 



n— s+p— 1 



+ 



+ P 1 S +Vn- S -l , 



(B.4) 
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and one gets the desired equation 



s + l 
p=0 



(B.5) 



because from (4.22) we have 



ps+l _ ps+2 _ i . (f + 1; n ) P s+1 _ ps+2 

^1 — ^1 — 1 ' ,,2 s+l ~~ s+2 

(s + l,n-(s + l)+p) ps+1 +1 _ +2 

..2 ^P ^P+l ~~ ^P+l • 

(B.6) 

The third of this relations is checked by splitting the last term of the sum at the 
end of P°Xl . With j p = s + 1 - ji - . . . - j p _! , it gives the factor (s+1 ' n ~jf +1)+p) 
multiplying an expression that can be rewritten as Pp +1 ■ The remaining terms in 
the sum are just PpXl ■ 

Now follows the proof of (4.24). Suppose that a = s + 1 > n + p — s = r + 1 = (3 
and write 



ps+l 

^p+l 



^ /T—p+l S—p+1 s 

-2J E (Ji,n + l-n)+ J2 (h.n+l-h)) 

^ \ j 1 = l j 1=r _ p+2 / 



X 



(»"-l-jl-...-jp-2 
E 0' 1 + -+ jp-l»W+p- 1 - jl - ...-jp-l) + 

S — l-jl— — -jp-2 v 

+ E ^ + - + Jp-i»w + p- 1 - ji - ... - j P ))x 

j P -i=r-ji-...-j p -2 
,r-ji-...-j p -i 

x( E 0'i + --+Jp^ + P-Ji---Jp) + 
V i„=i 

S-jl-----jp-2 

+ Yl (h + -+j P ,n + p-ji- ...-j p ) 

j'p=r+l-ji-...-j p _i 



(B.7) 
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The second term in the last factor is 



(r + 1, s) + (r + 2, s - 1) + . . . + (s - 1, r + 2) + (s, r + 1) = . 

(B.8) 

The biggest superior limit for the first term in this last factor, when multiplied 
by the second term in the previous factor, is determined by the value of j p -i = 
r + 1 — ji — ... — j p -2 , so we have that j p = and this product does not appear 
because ji = Vz. Following this argument, the same happens with all of the 
second terms, and only remains the product of the first ones that is nothing more 
than p£\ . 

Appendix C 

First we will proof (4.30). Taking into account (4.26) and (4.29), we have 

_ ..n{n— 1) 

P>0L,O>f3 

(C.l) 

n(n — 1 ) 

In what follows we omit the overall factor <=— ^ — . Because of (4.22), for a = (3 = n 
we get 

pn _ (ln)...(n-l,n) 



^2(n-l) 



(C.2) 



as required. 

For a < (3 = n, we have to prove that 

pn , f ^ ( a , n ) pn-1 , lf -,\n-(a+2) ( a , a +3)-( a , n ) pa+2 , 

r n+\ L ) ^ 2 + { L ) // 2(n-(a+2)) ^a+2 + 

, nw -(a+i) (a,Q + 2)...(a,n) pa+1 _ a (a, a + l)...(a, w - l)(a, n) pQ _ 

+ ^ ^ ^ 2 (n-(a+l)) ~+~ L > ^i 2 ^"") Q ~~ 

(C.3) 

In fact ^jr^-Pa = Pa+i because of (B.6), and the two last terms add up to 
, nn _( Q+ i) (a,a + 2)...(q,n- 1) (a + l,n) 

1 J jLt 2(n-(a+2)) ^2 > 

(C.4) 
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that can be expressed in terms of -P"+ 2 • Following this procedure step by step we 
arrive to 



(n-l,n) w _i_ 



(C.5) 

by (B.6). 

Thus, the n-th column of the kinetic matrix is null but for the element in the 
n-th row. Because of the symmetry that gives (4.24), we have for the n-th row 



pn , / .J (oi; n) „ n , , n _ (a+2) (a,a + 3)...(a,n) 

+ I i J ^ 2 + • • • + I // 2(n-(a+2)) ^a+2 + 

, , n n-(a+i) (a,a + 2)...(a,7i) pn , . n _ a (a, a + l)...(a,n) pn _ 

+ ^ ^ 2 (n-(a+l)) + I L ) „2(n-a) ^« ~~ U 



(C.6) 

where a < n indicates now the column. 

Take now f3 < a < n. Consider in (C.l) the terms with p = n, which are (C.6) 
with (3 in the place of a and give a null contribution. For a < p < n we can write the 
terms for a given p as 



K 



pa 



(P,n) 



2 {P,n-l){P,n) 



^ + ^-i(-l)^ + ^- 2 (-l) 

p p 

1)n -(0+i) (A^ + 2).-(g ,n) 
p" ' ' 

(p,p+l)...{p,n) 



+ ... 



p-n+(/3+l)^ ^ „2(n-03+l)) 



+ ^-n + ,(" 1 ) n " /3 



/i 



2(n-/3) 



(C.7) 



Applying (B.6), (C.7) becomes 
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K 



p p 

, rpp-i , (p- 1,^-2) pP -i u n2 {An-l)(^,n) 



n) 



/7/ 



, ! , (p -!,/?+!) ! , _ (/3+1) (/3,/3 + 2)...(An) 

• • • -l-l-^p-n+^+l) I" ^ 2 ^ p _ n+/3 |l ^(n-^+l)) 

,/pp-i , (p- pp-i lLl ^(MilHM' 



-7? (Z 9 - 1 ^) 
A* 



n) 



+ 



^P^" 1 / n n-(g+l) (A/^ + 2 )-(^^) , 
• -^p-n+^l L ) ^ 2 (n-(/3+l)) " + " 

pp-i , 



(C.8) 

The expression inside the brackets is the contribution to (C.l) from the terms 
with p — 1 divided by 1Z p -i j0 , . So, if the contribution for p is zero, it is zero too for 
p — 1 as long as (p — l,(3) ^ 0. When /3 < a, we have seen that for p = n we have a 
null contribution and then it is null for p = n — 1, and so is down to p = a + 1 and 
p = a. All the contributions being zero, the kinetic matrix is diagonal because of the 
symmetry. For (3 = a the same procedure works well until we reach p = a + 2, so that 
the contribution for p = a + 1 is zero, but here the chain stops because (p — 1, j3) = 
in the next step, and we get 



(n T iai) 



p 



n(n— 1) 



M 



E7"> pa 7? — 

'v aa -T (T _ n+a /Vera — 



p>Qf 



P 



n(n— 1) 



M 



■(-!)' 



(a, a + n) 

jj2(n-a) 



^+(-1)^^-1 + 
P 



2 (a,n-l)(a,n) pa 

+ l _i J ~4 r a-2 + 



P 



... + (-1) 



n — a+1 



(a, a + 2). ..(a, n) 



P 



2(n-(>+l)) 



pa 1 
a— n+a+1 1 



, _ a (a,a + l)...(a,n) 

V ' jj2(n-a) r a-n+a 
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(C.9) 



To get (4.30), we only need to prove that the square bracket in (C.9) gives 



(1, a). ..(a — 1, a) 

(CIO) 

To obtain this result, we apply to (C.9) the same procedure used in (C.8) and 



get 



(a — 1, a) 



pq-i |_ / ^ { a , n ) pq-i , l( ^ tl - a {a,a + l>...(a,w) Dtt _i 

r a.-\ ^\ L ) ^ 2 \ • • • \ \ J-J jj2(n-a) ^a-n+a-l 



(C.ll) 

We observe that a — n + a < a. If a — n + a < some of the terms in (C.9) are 
zero. If a — n + a = Z > we can repeat for (C.ll) the same operation we did on 
(C.9), and, after / + 1 steps, the last term is zero. After the step a — 2 we finally get 



(2, a). ..(a — 1, a) 

//2(a-2) 



Pi + (" 



-1)<^P? 



(1, a). ..(a — 1, a) 
//2(«-i) 



(C.12) 



as required. 



Now we show how to get (4.31). We will study the situation when (3 < a as long 



as the mass matrix is symmetric, and again omit 



Ml 



i(n-l) 



M 



, so we consider 



p>a,a>(3 



(C.13) 



Take a = n. Then, if j3 = n, one gets 



-MZ(l,n)(2,n)...{n-l,n) 

(C.14) 

for (C.13). If (3 < n it is the same expression obtained for the kinetic matrix with the 
factor — M% , so it vanishes. 

For a = n — 1 we have in (C.13) p = n, n — 1. The first value gives the same result 
as in the kinetic case but with the factor — M% , and it vanishes again for f3 < a. For 
p = n — 1 we have contributions with and without fx 2 : 
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pn / ^ (P,n) , D n / lN 2 <Aw-l)(^,n) 

^n-ll-J-J - ~2 h - r n-2l _i J "4 r 



+ ("^n-l) 



/// 



n— 3 



+ 



+ 



p n-i ( ir -g (A/3 + l)--(An) 



2(n-/3) 



(C.15) 

Inside the bracket for (i 2 , we have the same expression obtained for the kinetic 
part when p = n but without the first term, so it gives — P™/i 2 . Adding it to the next, 
we have 



(C.16) 

that, together with the rest of (C.15), is the same occurring in the kinetic case aside 
the factor — M 2 _ 1 . So it is zero if (3 < cm, and for f3 = a = n — 1 one has the result 
needed in (4.31). 

Similar rearrangements are used for the rest of the matrix. For example, for 
a = n — 2 one has two relations 

- K-lt + K-H-mD = Cr 2 2 (-M 2 _ 2 ) 

and one arrives again at the kinetic result with the factor — M 2 _ 2 . 
Appendix D 

Consider the four-derivative Lagrangian, 



(C.17) 



(D.l) 
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which also reads 



(D.2) 

where p = mfm 2 , and s = m\ + m 2 , . The covariant Ostrogradski method, in a 
slightly less refined version that the one presented in [9], would define a conjugate 
generalized momentum n = g^jjpj • The Legendre transformation performed on it 

leads to a Hamiltonian-like density from which the following two-derivative Helmholtz 
Lagrangian is derived 

C H = 7T ncf> + \ (12) tv 2 + 1 (12) 4> 2 + h JVC/) . 

(D.3) 

On the other hand, by using the auxiliary field technique of [17] the higher-derivative 
term is brought to second order by writing (D.2) as 

£4 = ~\ W) [p<p2 + S(p{n<p) + A(D0) " ' 

(D.4) 

where the equation of motion for A recovers (D.2) when substituted back in (D.4). 
Now, in spite of their quite different look, (D.3) and (D.4) are related by the simple 
field redefinition ix = — |-^y(A + scf>) . 
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